Draft version February 1, 2008 

Preprint typeset using rfTgX style emulateapj v. 10/09/06 



PROPERTIES OF A GAMMA RAY BURST HOST GALAXY AT Z 



P. A. Price, 2 A. Songaila; 



L. Cowie, 2 

KULKARNI, 6 



5 1 

D. B. 



Fox, 5 I. Hook; 



o 
o 

(N 



00 
Si 

6 



O 

in 
o 

o 



X 



J. Bell Burnell, 3 E. Berger, 4 A. Cucchiara, 5 
B. Penprase, 7 K. C. Roth, 8 B. Schmidt, 9 
Draft version February 1, 2008 

ABSTRACT 

We describe the properties of the host galaxy of the gamma-ray burst GRB060510B based on a 
spectrum of the burst afterglow obtained with the Gemini North 8m telescope. The galaxy lies at a 
redshift of z = 4.941 making it the fourth highest spectroscopically identified burst host. However, it is 
the second highest redshift galaxy for which the quality of the spectrum permits a detailed metallicity 
analysis. The neutral hydrogen column density has a logarithmic value of 21.0-21.2 cm -2 and the 
weak metal lines of Ni, S and Fe show that the metallicity is in excess of a tenth of solar which is far 
above the metallicities in damped Lyman alpha absorbers at high redshift. The tightest constraint 
is from the Fe lines which place [Fe/H] in excess of —0.8. We argue that the results suggest that 
metallicity bias could be a serious problem with inferring star formation from the GRB population 
and consider how future higher quality measurements could be used to resolve this question. 
Subject headings: cosmology: observations — galaxies: distances and redshifts — galaxies: abundances 
— galaxies: evolution — gamma-ray: bursts 



1. INTRODUCTION 

Long duration gamma-ray bursts are believed to form 
in the collapse of massive stars (e.g. Woosley 1993, 
Stanek et al. 2003) and hence to be a tracer of galaxies 
with ongoing star formation (Lamb and Reichert 2000). 
The afterglows associated with the burst are also ex- 
tremely bright immediately post-burst and can be seen 
to very large redshifts, currently out to z = 6.2 (Kawai 
et al. 2006). Prompt spectroscopy of the afterglow can 
therefore allow us to obtain the redshift of the galaxy in 
which the GRB has occurred and to study the metallicity 
of the interstellar medium in this host galaxy. 

With the caveat that we still do not fully understand 
the selection biases introduced by the formation of the 
GRBs, this may be the best way in which we can study 
the properties of individual galaxies at very high red- 
shifts (z > 5) in detail. The galaxies themselves are too 
faint (magnitudes of about 25) for detailed spectroscopic 
study, nor can they be identified as damped La absorbers 
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(DLAs) at these redshifts since the La forest becomes too 
thick to allow us to identify the DLAs and measure their 
column densities. 

However, in order to take advantage of GRBs to study 
the high redshift galaxy population we must be able both 
to localize the GRB quickly and then to rapidly identify 
it as a candidate high redshift object. The Swift satellite 
launched in 2004 (Gehrels et al. 2004) has made such 
studies possible by providing a large sample of GRBs 
with accurate positions to faint gamma ray detection 
thresholds. At the Swift limits Jakobson et al. (2006a) 
estimate that about 7% of the GRBs are at z > 5 so 
that there should be a handful of such objects available 
for study in each year of the Swift mission. 

Given this small number of available targets it is criti- 
cal to observe as many as possible with the highest qual- 
ity and highest resolution spectroscopy that we can ob- 
tain. However, even with the Swift data this remains 
extremely challenging since we must first identify the af- 
terglow and since the subsequent spectroscopy requires 
target of opportunity obervations (TOOs) on the largest 
ground based telescopes. Fortunately, some of the 8m 
class telescopes are operated in queue mode and can re- 
spond rapidly to such TOO events. The present program 
utilizes this capability on the two Gemini 8m telescopes 
where we can initiate observations with the GMOS spec- 
trographs (Hook et al. 2004) almost immediately when 
a candidate high redshift burst is identified. The GMOS 
spectral resolution (R~ 2000 in the mode which we used) 
is not as high as desirable but is adequate for measuring 
the redshift and providing metallicity estimates if we can 
obtain high signal-to-noise spectra. 

This Letter describes our GMOS observations of 
GRB060510B which we find to lie at a redshift of z = 
4.941. This is the fourth highest redshift identified for 
a GRB (e.g. Jakobson et al. (2006b) for a recent sum- 
mary) but is the second highest for which we can make 
a metallicity analysis. The highest redshift analysis is of 
the z = 6.2 GRB (Kawai et al. 2006). We find that the 
metallicity is in excess of a tenth solar, which is much 
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Fig. 1. — Spectrum of GRB060510B obtained with the 400 1/mm 
grating on the GMOS instrument on Gemini. We have made an 
approximately relative flux calibration of the spectrum based on 
the instrument throughput but have not attempted an absolute 
flux calibration. The exposure time was 4000s. The spectrum is 
shown in the rest frame of the host galaxy at a redshift z = 4.941. 
and the strong lines of La, Sill, OI, CII SilV and CIV are marked. 
The shaded regions show the position of the atmospheric bands 
and strong night-sky lines. 



higher than is seen in DLAs at z > 4 and also sug- 
gests that the metallicity of the GRB host galaxies is 
not changing as a function of redshift. 

2. OBSERVATIONS 

GRB060510B was detected by Swift on May 10th 2006 
at 8:22:24 UT. The XRT on Swift began observing the 
area 119s after the trigger and a bright X-ray source 
was located at RA=15 56 29.3 and Dec=78 34 09.4 
(J2000). The long duration gamma-ray burst was noted 
as a potential high redshift object (Krimm et al. 2006). 
A fading optical counterpart was detected by Mirabal 
and Halpern (2006) at RA=15 56 29.615 and Dec=78 
34 13.02. The red (R-I) color was again suggestive of 
a possible high redshift candidate. NIRI observations 
with Gemini-North detected the counterpart with J~19 
at just over 2 hours after the trigger (Price et al. 2006a). 

We obtained a spectrum with the GMOS spectro- 
graph on Gemini North, commencing at 2006 May 10.547 
UTC just over 2.5 hours after the burst (Price et al. 
2006b). Four 1000s exposures were obtained with the 400 
lines/mm grating giving a resolution of 1900 and wave- 
length coverage from 5950 A to 10200 A. Standard CCD 
reduction steps were performed with the gmos package 
within IRAF, pairs of exposures were sky-subtracted and 
combined, the spectra extracted using apextract, and 
the two spectra summed to yield the final product. The 
spectrum was approximately relatively flux calibrated us- 
ing the instrument throughput. The final spectrum is 
shown in Figure [TJ 

The spectrum shows a damped La line and an exensive 
set of metal lines at a redshift of z — 4.941 which we 
take to be the redshift of the host galaxy. Some of the 
stronger features in the spectrum are marked on Figure[TJ 
The spectrum also shows a strong break across the La 
wavelength corresponding to the onset of the La forest 
at shorter wavelengths. The mean transmission in the 
La forest is 18%, consistent with La forest transmissions 
seen in quasar spectra at this redshift (Songaila 2004). 
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Fig. 2.— Damped La and L/3 profiles for N(H1) of 1 x 10 21 
cm -2 , 2 X 10 21 cm -2 and 4 X 10 21 cm" 2 . 



3. HOST PROPERTIES 

We first measured a redshift of z — 4.941 from the 
weaker singly ionized metal lines in the spectrum. Fits to 
the damped La and L/3 lines centered at the wavelength 
corresponding to this redshift are shown in FigureO The 
best fit to the red wing of the damped La profile is given 
for a N(HI)=1.7x 10 21 cm -2 and higher column densities 
are prohibited. The L/3 profile favors a slightly lower 
value of N(HI) = 1.0 x 10 21 cm -2 but this constraint is 
based on a single pixel which could be contaminated. 
We therefore adopt a range of log/V(Hl) = 21.0 - 21.2 
for the logarithmic column density. 

A second strong absorption feature seen in the La 
forest at 7100 A is not saturated at L/3 and is not a 
DLA. This neutral hydrogen excess may be caused by the 
higher density intergalactic medium in the vicinity of the 
galaxy and it will be interesting to see if this is common 
in GRBs at these redshifts since it could provide consid- 
erable infomation on the structure and ionization of the 
overdense regions in which the galaxies are forming. 

Although we detect many lines from the galaxy, most 
are strong and, in low resolution spectra like the present 
one, the lower limits on the column densities which can 
be obtained from such lines are not particularly useful. 
We therefore focus on the weaker lines of Ni, S and Fe 
which are seen in the spectrum. The absorption profiles 
of some of these lines are shown in Figure [H 

Even in these moderately weak lines saturation can be 
a problem and saturated components can be masked by 
unsaturated features in the lines. This problem, long fa- 
miliar in interstellar medium studies, has recently been 
re-emphasized by Prochaska (2006) where an extensive 



z ~ 5 GRB host galaxy 



3 




4.9410 

1 ■ 1 



-2000 -1000 1000 2000 
VELOCITY (KM/S) 




4.9410 

-2000 -1000 1000 2000 
VELOCITY (km s" 1 ) 

Fig. 3. — Upper panel, bottom to top: lines of SII1251, SII1254, 
FeII1608, FeII1260. Lower panel, bottom to top: lines of N5II1317, 
NiII1370, NiII1455 and NiII1467. For each species the lines are 
shown in order of decreasing oscillator strength with the stronger 
lines at the bottom and the weaker at the top. The higher noise 
at the Fel608 line is a result of decreasing throughput at longer 
wavelengths. 

discussion and critique of GRB host metallicity measure- 
ments can be found, as well as the many historical ref- 
erences on the topic. For the present work we measured 
the various Ni, S and Fe lines under the assumption that 
they are unsaturated and lie in the strictly linear portion 
of the curve of growth. The column densities measured 
in this way and their one sigma uncertainties are given 
in Table 1. We emphasize that these are only lower lim- 
its to the column density. If there are narrow saturated 
components concealed in the absorption lines the column 
density could be higher. However, even these lower limits 
already place a lower limit on the meta llicity of about — 1 
relative to solar (Asplund et al. (2005) ; Table 2) which 
is already considerably higher than the values seen in 
DLAs at z > 4 (Prochaska 2003). We illustrate this in 
Figure [H 

If the signal to noise is high enough we can use very 
weak lines to measure upper bounds on the metallicity 
even in low resolution spectra (e.g. Savaglio 2006) . How- 
ever, the quality of the very high redshift GRB spectra 
is not generally adequate to do this. We tested the limits 
we could obtain in the present data using the very weak 
Fel611 and Nil467 lines. In each case we created a model 
with two components each with b — 8 km s _1 and two 
broader (b — 50 km s" 1 ) components and adjusted the 
column densities to fit the stronger lines. (Here b is 0.60 
times the full width at h alf m aximum of the absorption 
line (Cowie & Songaila (|l986h ). We varied the column 
densities in the narrow components and compared the 



TABLE 1 

Column Densities 



Ion 


A (A) 


/ 


log N (cm 2 ) 


HI 


1215.67 


0.416 


21.1 ±0.1 


Fe II 


1608.45 


0.062 


> 15.9 


S II 


1250.58 


0.005 


15.6 (14.8 - 15.9) 




1253.81 


0.011 


15.5 ±0.2 


Ni II 


1317.22 


0.145 


14.4 ±0.2 




1370.13 


0.131 


14.0 ± 0.2 




1454.84 


0.0595 


14.4 (14.0 - 14.6) 


Note. — 


All errors 


are la. 








TABLE 2 








Metallicities 




Element 




Solar 


Relative 


S 




-4.86 


-0.8 ±0.3 


Fe 




-4.55 


-0.3 to -0.8 


Xi 




-5.77 


-1.2 ±0.3 



resulting Voigt profiles with the weak lines. Within the 
noise and continuum fitting uncertainties we cannot ob- 
tain a robust upper limit on the column densities. 

The spectrum also shows strong high ionization lines 
with logiV(NV) in excess of 14.5 and logiV(CIV) and 
logiV(SilV) in excess of 15. Strong fine-structure lines 
are seen for Si II, C II and O I, as i s com mon to all 
the GRB absorbers (Vr eeswii k et al. (|2004h : Berger et 
al. (|2QQ6T ): Chen et al. (|2005[ )). Given the uncertainty 
in estimating the ground state column densities we do 
not attempt to infer gas densities but it is clear, as has 
been noted by many authors, that GRB absorption is 
produced in much higher gas density enviroments than 
DLAs in quasar sightlines, which only extremely rarely 
show such lines (e.g. Chen et al. l|2005h ). 

4. DISCUSSION 

In Figured] we show all of the currently available metal- 
licity measurements for z > 2 GRB host galaxies from 
the references summarized in the figure caption. The 
present measurement of GRB060510B is shown as the 
filled diamond with upward pointing arrow to emphasize 
that it is strictly a lower limit. We have shown other 
measurements of the GRB hosts as open diamonds. We 
show measurements based on low resolution and moder- 
ate signal to noise spectra with upward pointing arrows 
on 1 <T error bars while values based on high resolution 
or high signal to noise observations are shown with just 
the 1 a errors. Irrespective of redshift the bulk of the 
GRB metallicities reported so far lie at or above about 
a tenth solar. However, GRB hosts clearly have a wide 
range of metallicities, for example, compare GRB050730 
(z = 3.969; Chen et al. 2005) with a metallicity of -2 
and GRB000926 (z = 2.038; Savaglio et al. 2003) with a 
value of —0.13. 

The typical metallicity of GRB host galaxies is con- 
siderably higher than the metallicities found in DLAs, 
which are shown in Figured] as the filled squares. How- 
ever, there is only one z > 4.5 DLA with a measured 
metallicity (Songaila and Cowie 2002) so the comparison 
sample at the highest redshifts is very limited. This ef- 



4 



Price et al. 



x 

LJJ 




2.0 2.5 3.0 



3.5 4.0 4.5 5.0 
REDSHIFT 



5.5 6.0 6.5 



Fig. 4. — Metallicity relative to solar versus redshift for the GRB 
host galaxies (diamonds) superimposed on the corresponding mea- 
surements of metallicity in damped Lyman alpha systems along 
quasar sightlines, taken from the compilation of Prochaska et al. 
(2003). The value measured in this paper is shown as the filled dia- 
mond and is based on the SIT, Fell, and Nill measurements of Table 
2. The open diamonds show values from the literature (Savaglio et 
al. 2003, Vreeswijk et al, 2004,2006 Chen et al. 2005, Starling et 
al. 2005, Watson et al. 2005, Berger et al. 2006 and Fynbo et al. 
2006). For lower resolution measurements (including the present 
measurement) where saturation effects could significantly raise the 
measured metallicity we plot the data with an upward pointing 
arrow. The dashed line shows the metallicity above which GRBS 
are not expected to form in collapsar models. 

feet and the higher gas densities evidenced by the fine 
structure lines are natural consequences of the selection 
biases. GRB sightlines target star forming regions of 
galaxies where densities and metallicities will be higher, 
whereas the cross section- weighted DLAs in quasar sight- 
lines probe more extended lower density and lower metal- 
licity regions in galaxies and may also be weighted to 
intrinsically lower luminosity galaxies. 

Ultimately we would like to use GRBs to probe the star 
formation history of the universe (Lamb and Reichert 
2000, Price et al. 2006). They have many advantages for 
this purpose not the least of which is that they sample 
the rate at which individual stars form irrespective of the 
mass or luminosity of the host galaxy. We can therefore 
measure the total star formation over the entire range 
of galaxies. However, before we can do this we must 
understand the biases in the selection of the stars which 
become GRBs and determine whether we can succesfully 
allow for such selection effects. 

The most probable bias is that GRBs occur only in low 
metallicity galaxies and that we may therefore miss all 



of the star formation occuring in more evolved systems. 
In the collapsar models GRBs are formed in single mas- 
sive stars only if the metallicity is below ~0.3 Zq (e.g. 
Woosley and Heger 2006). If such a bias is present it 
could introduce a strong redshift dependence in the in- 
ferred star formation history if the average metallicity of 
the host galaxies is lower at high redshift. 

The evidence of Figure [4] is somewhat confusing as re- 
gards this issue. In the first place it shows that while 
most metallicities in GRB hosts could fall below 0.3 
Z©, at least one (GRB 000926, Z = 0.7 Z , Savaglio 
et al. l|2003h ) has a metallicity well above the critical 
value. Fynbo et al. (2006) suggest that this indicates 
that collapsars resulting from single massive stars are 
not the only progenitors of long GRBs or that massive 
stars with Z > 0.3 Z Q can also produce long GRBs. It is 
also possible that this particular GRB formed in a lower 
metallicity region of the galaxy than that traversed by 
the sight-line. However, if some of the measured metal- 
licities are underestimates this problem would become 
much more severe and we would have to conclude that 
GRBs regularly formed in high metallicity environments. 
This issue clearly requires a much larger sample of very 
high quality spectra to resolve. 

Furthermore if we do assume that GRBs are primarily 
formed in lower metallicity stars then Figure |4] strongly 
suggests that metallicity bias will be extremely impor- 
tant in determining star formation rates. Most of the 
measured values lie just below the critical value and pre- 
sumably as their metallicity continues to increase would 
drop from the sample. Since a large fraction of the star 
formation would occur in the higher metallicity objects 
we would miss much of the activity. With precise metal- 
licity measurements we could make a detailed compari- 
son with models of the star formation to see if the distri- 
bution was consistent and the metallicity threshold was 
evident but such a comparison is premature given the 
present quality of the data. 
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